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SUMMARY 


Formulas are presented for the determination of the creep deflection- 
time characteristics of an initially curved idealized H-section column. 
These results were ottained from closed-form solutions of the differen- 
tial equation of tending (derived in KACA TN 3137) of a beam column whose 
creep properties .are of a nonlinearly viscoelastic nature. The critical 
time (the time required for infinite deflections to develop) established 
by these solutions is tabulated and plotted for a wide range of the 
parameters involved. 


INTRODUCTION 


The effect of creep on the behavior of initially curved columns was 
previously investigated in reference 1, and the resulting deflection- 
time characteristics were obtained for several values of the parameters 
involved. It was found that every column whose material is subject to 
nonlinear creep — and this includes all columns made of structural metals 
such as aluminum, steel, titanium, and so forth, when subjected to high 
temperatures — buckles if the axially compressive load acts upon it for a 
sufficiently long time. This statement is true even though the compres- 
sive force is less than the static critical load which is defined as that 
load which would cause buckling instantaneously. Depending upon the ratio 
of the applied load to the static critical load, the initial deviation 
from straightness, and the creep properties of the metal, the time 
required for the development of infinite deflections — the so-called 
critical time — may be anywhere between a few seconds and a few years. 

It is of great importance to the structural designer of supersonic air- 
craft to know how much the critical time of his structure is. 

In the present report the differential equations derived in refer- 
ence 1 for the analysis of the behavior of idealized H-section columns 
are solved in closed form for integral values of the exponent in the power 
function defining the assumed creep law. These solutions are used for the 
determination of the deflection as a function of time of an initially 
curved column whose end load is less than the static critical load of the 
col umn . The critical time is calculated for a wide range of the exponent 
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as a function of a parameter which includes the effect of end load and 
initial curvature. The results of the calculations are presented in 
tables and charts which enable the designer to determine the critical 
time once the parameters appearing in the basic uniaxial tensile or com- 
pressive creep law are determined experimentally at design temperature. 

The authors are indebted to Professor N. J. Hoff for his guidance 
and criticism, to Professor F. V. Pohle for his assistance on the mathe- 
matics contained herein, to Mr. R. Ortasse for his assistance with the 
calculations, and to the National Advisory Committee for Aeronautics for 
sponsorship of the research reported in this paper. 
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area of idealized H-section 
constants 

effective elastic modulus 

amplitude of time -dependent deflection (accrued for 
t > 0) divided by h 

amplitude of initial deviation from straightness of 
unloaded column divided by h 

amplitude at t = 0 of total deviation from x-axis of 
loaded column divided by h (x-axis is drawn between 
end points of column) , f^/^l - crg^j 

distance between flanges of idealized H-section 

moment of inertia of idealized H-section, Ah 2 /t 

integrals 

column length ■ 

exponent in viscosity term 

integer 

integer 


Cm-D/2 

S = (-l) n cos m_2 (nir/m) 

n=l 
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t time 

t cr critical time 


x axial coordinate 

z amplitude for t s? 0 of total deviation from x-axis 

divided by h, f c + frp^ 


€ 

*o 

X 

a 

a 

ct E 

°o 

r 



$1* $2 


strain 

strain at t = 0, a 0 /E-]_ 

viscosity coefficient 
stress 

average axial compressive stress 
static truckling stress, Jt^E-^l/AL^ 


constant stress 


time parameter. 




critical value of time parameter ( t » r cr as 

time parameter corresponding to z = l/2 
functions of x 


z 


(') ■*( )/8t 
( )' )/Sx 


00) 


CHEEP LAW 


In the analysis of reference 1 the fundamental uniaxial tensile or 
compressive constant- stress creep curve is approximated by a straight 
line, the slope of which can be considered as the secondary creep rate of 
a real material (fig. l) . This idealization accounts for the actual second' 
ary stage of creep (assumed to be of a viscous nature), approximates the 
initial elastic or elastoplastic stage and the primary creep stage, and 
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ignores any final stage. The relation between the strain rate and the 
stress corresponding to the simplified creep curve is 

6 = (ff/Bj + (t^A) (1) 


in which e and a, respectively, are the strain and stress, Ej is the 

effective modulus (fig. l), m and X are parameters defining the vis- 
cous behavior of the material, and the dot over a symbol indicates dif- 
ferentiation with respect to time t. The three material parameters E]_, 

m, and X are considered as constants for a given temperature and can be 
determined experimentally from conventional tensile or compressive creep 
tests. The corresponding relationship among stress, strain, and time for 
such tests can be found from equation (l) if the strain g Q at t = 0 

is taken as cr 0 /Ei (fig. l). Thus 



(E]/X)c 


m-1 


+ 1 


( 2 ) 


DEFLECTIONS OF A COLUMN WITH INITIAL CURVATURE 


Differential equations were derived in reference 1 defining the 
deflection-time characteristics of a simply supported idealized H-section 
column (figs. 2 and 3) whose material parameters are the same for tension 
and compression. The column is assumed to be loaded instantaneously at 
t = 0 with an axial load which remains constant for t >0. The differ- 
ential equations derived were simplified with the aid of the assumption 
that the sinusoidal shape of the centroidal axis of the unloaded column 
was maintained in the loaded column and only the amplitude varied. The 
differential equations so obtained were readily solved in integral form. 
Solutions to these equations are presented in table 1 in which the equa- 
tions referred to are from reference 1. 

In appendix A the integrals appearing in table 1 are evaluated In 
closed form for any even or odd integral value of m. Hence the following 
relations govern the deflection-time characteristics of the column ana- 
lyzed for a < cr E . For m an odd integer (m > l) with 0 < fm < 

■^o 

and frp ^ z ^ 00 , 

±0 - ~ 


00 
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- ( z: /a (-xf 

n=l 


cos ) 


log + W^J • 

> 


For m an even integer (m >0) with 0 < fm ^ l/2 and frj ^ z ^ l/2. 


■ ^A 108 (t) 


( m/ 2 ) — 1 


n»l 


(-l) n cos m_2 ^) log 


j\z 2 + tan 2 ^ 

f3 , 

[4f To 2 + tan 2 

(¥)] 


(5) 


/ 


T, = 


( 2 ^) (log /_i_\ + ^x:' 1 (-D n cos m " 2 (^) log- [ x ^ tan (mjl 

' m \ 2f T 0 / 11=1 [^ f T 0 + tan (¥ 


m] 

C 6) 


with 0 < fip ^ l/2 and l/2 < z ^ eo, 


-1 * (*£) n £ *“®3 - @} 
n-X, J ; • • • L J 


C7) 
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T cr = t-]_ + ( 2 m “ 2 Tt/m 2 ) y ( -l)^ m-n-1 ^/ 2 (m-n) sln m " 2 (nir/2m) (8) 

n=lj 3 j • ■ • 


with l/2 fjp Q < 00 and frp^ ^ z < co, 



( 10 ) 

If m = 1, the corresponding column is constructed of a linearly 
viscous material (Maxwell material) and from reference 1 its exact 
deflect ion- time behavior is governed by the relations (see also refs. 2, 
3, and b) 


z 


fm e 
■“-o 


2T 


( 11 ) 


and 


'cr 


(11a) 


Equations (3) to (lla) define completely the deflection-time charac- 
teristics of col umn s whose properties have been previously defined. It 
may be noted that for all values of m greater than unity there exists a 
finite critical time which is defined as that time at which the column 
deflections increase. without limit. 
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DISCUSSION 


The results presented in equations (3) to (lla) can "be used to deter- 
mine the deflection as a function of time of col umn s whose properties are 
approximated by those described earlier. With the aid of these equations 
curves of deflection versus time for various values of the zero-time 
deflection parameter f^ can be computed. Such curves are presented 

in reference 1 for a wide range of values of fij and for m = 1, 2, 3 , 4, 

and 5- 

Since the critical time parameter T cr is a measure of the life 

span of a column, this quantity is perhaps the most significant parameter 
of the present analysis. Hence in table 2 and figures 4 and 5 T cr is 

given as a function of fij Q for a wide range of the exponent m. The 

results for integral values of m were obtained with the aid of equa- 
tions (4), (8), and (10). However, for large values of m it is more 
convenient to determine r cr from numerical integration of the pertinent 

relations given in table 1, than it is to apply the corresponding closed- 
form solutions. It may be noted that results for m = 1.5 and m = 1.1 
are also included in table 2 and figure 4. These latter results were 
obtained by performing the integrations indicated in table l(a) in closed 
form for m = 1.5 and numerically for m = 1.1. In order to obtain the 
actual critical time t cr of a given column with known values of the 

applied stress a and amplitude of initial deviation from straight- 
ness f^, results of uniaxial tensile or compressive creep tests corre- 
sponding to the design temperature must be available. From such infor- 
mation the parameters Ej_, m, and X can be determined. If these 

parameters differ significantly for tension and compression tests, it is 
suggested that they be chosen to correspond to the compression tests, 
since for all total -deflect ion amplitudes zh < h/2 both flanges of the 
H-section column are in compression. In view of the simplifying assump- 
tions regarding the creep law, shape- of cross section, and affineness of 
the shapes of the unloaded and loaded column,, the present results can be 
considered only as a first approximation in the analysis of actual columns. 


Polytechnic Institute of Brooklyn, 
Brooklyn/ N. Y., Dec. 2, 1952. 
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APPENDIX A 

EVALUATION OF INTEGRALS APPEARING IN 
DEFLECTION-TIME RELATIONSHIPS FOR 
INTEGRAL VALUES OF EXPONENT m 

The integrals required for the evaluation of the various deflection- 
time relationships arising from the analysis of the creep deflections of 
an axially loaded initially curved bar are from table 1 

1 1 = fdz/fcl/2) + z] m - [(l/2) - z] m ’ (Al) 

and 

1 2 C z + + [ z " (l/2)] m | (A2) 

For integral values of m, it is seen from table 1 that 12 requires 
evaluation for both odd and even values of m, whereas I 2 need be 
determined only for even values of m. 

In the ensuing calculations, the integrals Iq and I 2 are replaced 
hy . . 

Jq = jfix/fcqCx) _ (A3).. 

and 

J 2 = j/dbc/a> 2 (x) (A^) 

in which ®q(x) = j(l + x) m - (l - x) 111 ], ® 2 ( x ) = [(x + l) m + (x - l)^, 

x = 2z, Jq = 2^ _m Iq, and J 2 = 2^ -m I 2 . Each of the required integrations 
will be performed with the method of partial fractions. 
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Determination of Iq for Odd Values of m 

If the denominator of the integral in equation (A3) is equated to 
zero, then the following relationship is obtained for the roots x of 
the resulting odd mth-order polynomial: 


(1 + x)/(l - x) = l 1 ^ 


(A5) 


in which the required roots of unity are 

l l/m - e< 2 " r /“> 1 , n = 0, ±1, ±2, . . ±(».- l)/2 (AS) 


Hence, substitution for 


x l/m 


from equation (A 6) into equation (A5) and 


subsequent solution- for x yield 

x = i tan (nir/m), n = 0, ±1, +2, . . . , ±(m - l)/2 (A7) 

Thus, the denominator of equation (A3) can be factored as follows: 

$q(x) = xjx 2 + tan^^jjj^x 2 + tan 2 ^^JJ . . . . x 2 + tan 2 j^^^(m - l)/2 j* 

(A8) 

The integrand of Jq can be resolved into partial fractions to yield 


l/ <I>q(x) = (A/x) + 


U-l)/2 


n=l 


(B n x + C. 


n)/E^ + tan 2 (njr/m)J (A9) 


in which A, B n , and Cq are constants. From equation (A9) 


®l(x) ^ 


(m-l)/2 . 

(A/x) + ( B n x + ^n)/L x2 + tan 2 (nrt/m)J 

n=l ' 


^ = 1 (A10) 


and hence A is determined from the condition that, since $q(o) = 0 


lim $q(x)( A/x) = lim A$q'(x) = 1 
x-*0 x-»0 


(All) 
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in which the prime indicates differentiation with respect to x. 


A = l/2m 


Similarly, since tan (nir/m)^ = 0, 


(m-l)/2 

lim ®l(x) (®n x + ^n)/!-? 2 + tan 2 (nTt/m)J 

x->i tan (nir/m) n=l 


= lim 4 1 ’(x)(B n x + C n )/2x 

x->i tan (nir/m) 


= 1 


for each n = 1> 2, . . . , (m - l)/2. Thus 

|jBni tan (nir/m) + C^]/2i tan (njr/m)| ^'[i tan 
From e quat ion ( A7 ) 

(1 ± x) = e^H^^/cos (ntf/m) 


= 1 


and hence 

®^'[i tan (n«/m^3 = 2m(-l) n cos 2-m (n3t/m) 
Equations (Al4) and (Al6) yield 

^ = (l/m)(-l) n cos m-2 (nrt/m) 

and 


Hence 

(A12) 

(A13) 

(Al4) 

(A15) 

(Al6) 

(A17) 


Cn = 0 
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From equations (A 9), (A12), and (A17), the integrand of J x becomes 


, (m-l)/2 „ _ o 

l/®l ( x ) = (l/2m)(l/x) + (l/m) (-1) cos (nir/i 

n=l 


m) jx/js 


x + tan^nir/i 


4 


(A18) 


Hence from equations (Al), (A3), and (Al8) , after performance of the 
indicated integration, the following expression is obtained for the indef- 
inite integral I]_ for odd integral values of m > 1 


I x = ( 2 m " 2 /m) - 

(m-l )/2 n m 2 r 1 

log 2z + 'y (-l) cos (mr/m) log |_4z^ .+ tan^(n 7 t/m)j * 


n=l 

^ J 


(A19) 


Evaluated between the limits z and fm , I, becomes 

± 0 


I 


1 



Also 



(A20) 



lim Ii(z,f T ) 

z— »« 


( 2 "~ 2 V 

log 

( 1 \ 

- (u z: /2 (-ir =o S - 2 (^) io g [ 

n=l L 

tf T 0 2 + tan2 (ii 

til 

\ m ) 

c 

CVJ 



s. 






(A21) 


In the determination of equation (A2l) use was made of the relationship 


(-l) n cos m ^(mr/m) = -l/2 

n=l 


( A22) 
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in which m is an odd integer greater than unity. This relation is 
derived in appendix B. 7' 


Determination of Ip for Even Values of m 

The present analysis is performed in a manner analogous to that of 
the preceding calculations. Since m is now an even integer the poly- 
nomial <B,(x) in equation (A3) is of order m - 1. The required 

(m - l) roots of <5p(x) are then 

x = i tan (nff/m), n = 0, ±1, ±2, . . ±[(m/2) - £} (A23) 

Since equations (A7) and (A23) differ only in the choice of values of n, 
the resulting expression for. Ip of the preceding section can be used 

in the present calculations provided the upper limit on the summation 
sign in equation (A19) is replaced by (m/2) - 1. Hence for even integral 
values of m 
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Since 


x i 1 - -ie ±f ” /2B)i /= ln C«/am) 


(A31) 
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l4 


and 

®2'L"^ co ^ (x«/ 2 n)] = 2mi^ m (-l)^ n sin 2 m (mt/2m) (A32) 

B n = 0 

and 

C n = (l/m)(-l)^ in n cos (mt/2m) sin™ ^(mt/2m) (A 33 ) 

Therefore .from equations (A2), ( A4) , (A30), and (A33) for even integral 
values of m 


in X 

12 = ( 2 m- 1 /m) > ( _]_)( m-n-1 )/ 2 sin m_ 2 (n 7 t/ 2 m) tan tan (nn/2m)\ 

n=l,3 ) 5> • • • 

(A34) 




(A3T) 
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I 2 (“,l/2) = (2 m_2 Vm 2 ) XZ (-l) (m " n_l)y/2 (m-n) sin m " 2 (n^/2m) 

n =l >3,5, • • • 


(A38) 
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APPENDIX B 


( m-l) /2 p 

DERIVATION OF THE RELATION (-l) cos^^mt/m) = -l/2 

n=l 

FOR ODD INTEGRAL VAIUES OF m 


In appendix A I-^Mjftp^ was determined with the aid of the 
t ions hip 

(m-l)/2 m-2. 

S = > (-1) cos (nrr/m) 

n=l 


= - 1/2 

in which m is any odd integer greater than unity. Equation (Bl) 
now be shown to be valid. From reference 5,, page 69, 

. Cm-3)/2, m 0 . 

co6 m-2 (nrt/m) = 2^ -m y ( r ) cos C m " 2r - 2)(nrt/m)J 

r=0 V 

in which the binomial coefficient (ref. 5, p. 19) is defined as 


/ m-2\ _ / m-2 \ _ 

\ r ) ~ \m-2-rJ r . 


(m - 2)1 


!(m - 2 - r) .’ 

Equations (Bl) and (B2) yield 
(m-l)/2 (m-3)/2 

S = 23" m 

n=l r=0 

With the order of summations interchanged, equation (b 4) becomes 
{“ill/ 2 ,- I/ 2 . ,n 


l 

(-l) n ( m r ^) cos [( m - 2r - 2)(nn/m)] 


S = Z 


3-m 


( m r 2 ) ZZ (~l) n cos [(m - 2r - 2)(jt/m)n] 
n=l 


r=0 


rela- 


(Bl) 

will 

(B2) 


(B3) 

W) 

(B5) 
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From reference 5, page 82, 
(m-l)/2 


m-l)/2 

(-1) cos I ( m - 2r - 2)(ir/m)n] = 
n=l L 

-(l/2)|l - (-l)^ m_1 ^ 2 cos £(m - 2r - 2)( it/2)J/cos Qm - 2r - 2)(Jr/2m)jj , 

(B6) 

Since m is odd, m - 2r - 2 is also odd, and therefore equation (b 6) 
reduces to 


(m-l)/2 


n=l 


(-l) n cos Jjm - 2r — 2}(ir/m)nJ = -l/2 


(B7) 


Equation (B7) together with equation (B5) yields 


S = -2 


2-m 


(m-3)/2 


r=0 


(”; 2 ) 


(B8) 


From reference 5, page 19, 


m-2 




= 2 


m-2 


(B9) 


The number of terms under the summation sign in equation (B9) is m - 1 
which, since m is odd, is an even number. From equations (B3) it can 
he seen that the binomial coefficients are symmetric. Hence the sum of 
the first (m - l)/2 terms under the summation sign of equation (B9) is 

(l/2)(2 m_2 ) = 2 m_ 3. in order to obtain this sum from equation (B9), let 
r = 0, 1, . . ., (m - 3)/2. The number of terms considered is then 
(l/2)(m - 3) + 1 = (m - l)/2. Consequently 


(m-3)/2 


r=0 


rf) - 


,m-3 


(BIO) 


and, therefore, equation (B8) yields 

S = -1/2 


which was to be proved. 
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TABLE 1.- INTEGRALS REQUIRED FOR DETEBMXNATIOR OF t AS A FUNCTION OF 
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[Reproduced from ref. 
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(33). 
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TABLE 2.- VALUES OF T cr FOR 1.1 £ m £ l4 AND 0.01 ^ f T % 1.00 



1.1 

1.5 

2 

3 

4 

5 

6 

7 

0.01 

7.165 

3.178 

2.741 


3.891 

5.612 

8.502 

13.73 

.02 

6.918 

2.852 

2.395 

1 

3.198 

4.505 

6.759 

10.58 

.05 

6.471 

2.420 

1.937 

■SSI 

2.286 

3.053 ' 
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.0940 

.70 
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.0557 
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.144 
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.0344 
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.0633 
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.0220 

1.00 

4.997 
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,464 

.187 

.0911 

.0477 

.0259 

.0145 
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8 

9 

10 

11 

12 

13 

l4 
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197-4 
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623.6 

wmm 
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.07 
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17.79 
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71.36 
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.10 
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8.018 

11.07 

16.59 

25.ll 

38.33 

58.95 

.20 

1.716 

2.169 

2.754 

3.560 

4.598 

6.021 

7-940 

.30 

.667 

.745 

.828 

• 932 

1.058 

1.213 

1.400 

.4o 

.292 

.291 

.287 

.288 

.290 

• 295 

.304 

.50 

.143 
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.111 

.100 

.0909 

.0833 

•.0769 

.60 

.0733 

.0583 

.0469 

.0386 

.0319 

.0265 

.0223 

.70 

.0398 

.0291 

.0215 

.0161 

.0122 

.0094 

.0072 

.80 

.0227 

.0153 

.0105 

.0073 

.0051 

.0036 

.0025 
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.0134 

.0085 

.0053 

.0035 
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Figure 1.- Idealized creep curve. 
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Figure 2.- Idealized H-Bection. 
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Figure J.- Deflections of simply 
supported column. 
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